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Abstract — The design of modulation schemes for the physical 
layer network-coded two way relaying scenario is considered with 
the protocol which employs two phases: Multiple access (MA) 
Phase and Broadcast (BC) Phase. It was observed by Koike- 
Akino et al. that adaptively changing the network coding map 
used at the relay according to the channel conditions greatly 
reduces the impact of multiple access interference which occurs 
at the relay during the MA Phase and all these network coding 
maps should satisfy a requirement called the exclusive law. We 
highlight the issues associated with the scheme proposed by 
Koike-Akino et al. and propose a scheme which solves these 
issues. We show that every network coding map that satisfies the 
exclusive law is representable by a Latin Square and conversely, 
and this relationship can be used to get the network coding maps 
satisfying the exclusive law. Using the structural properties of 
the Latin Squares for a given set of parameters, the problem 
of finding all the required maps is reduced to finding a small 
set of maps for A/ PSK constellations. This is achieved using 
the notions of isotopic and transposed Latin Squares. Even 
though, the completability of partially filled M x M Latin 
Square using A/ symbols is an open problem, two specific cases 
where such a completion is always possible are identified and 
explicit construction procedures are provided. The Latin Squares 
constructed using the first procedure, helps towards reducing the 
total number of network coding maps used. The second procedure 
helps in the construction of certain Latin Squares for A/-PSK 
signal set from the Latin squares obtained for A//2-PSK signal 
set. 

I. Background and Preliminaries 

We consider the two-way wireless relaying scenario shown 
in Fig[T] where bi-directional data transfer takes place between 
the nodes A and B with the help of the relay R. It is 
assumed that all the three nodes operate in half-duplex mode, 
i.e., they cannot transmit and receive simultaneously in the 
same frequency band. The relaying protocol consists of the 
following two phases: the multiple access (MA) phase, during 
which A and B simultaneously transmit to R and the broadcast 
(BC) phase during which R transmits to A and B. Network 
coding is employed at R in such a way that A (B) can decode 
the message of B (A), given that A (B) knows its own message. 

The concept of physical layer network coding has attracted a 
lot of attention in recent times. The idea of physical layer net- 
work coding for the two way relay channel was first introduced 
in |[T|, where the multiple access interference occurring at the 
relay was exploited so that the communication between the 
end nodes can be done using a two stage protocol. Information 
theoretic studies for the physical layer network coding scenario 
were reported in f2^, f3|. The design principles governing the 
choice of modulation schemes to be used at the nodes for 
uncoded transmission were studied in Q. An extension for 
the case when the nodes use convolutional codes was done 
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Fig. 1. The Two Way Relay Channel 



in fSl. A multi -level coding scheme for the two-way relaying 
scenario was proposed in [6|. 

It was observed in Q that for uncoded transmission, the 
network coding map used at the relay needs to be changed 
adaptively according to the channel fade coefficient, in order 
to minimize the impact of the multiple access interference. 

A. Signal Model 

Multiple Access (MA) Phase: Let S denote the symmetric 
M-PSK constellation used at A and B, where M = 2^, X 
being a positive integer Assume that A (B) wants to transmit 
an A-bit binary tuple to B (A). Let : S denote the 

mapping from bits to complex symbols used at A and B. Let 
— 1^{sa), xb = G ^ denote the complex symbols 

transmitted by A and B respectively, where sa,sb G ^2'"- The 
received signal at i? is given by, 

Yr = Haxa + Hbxb + Zr, 

where Ha and Hb cue the fading coefficients associated with 
the A-R and B-R links respectively. The additive noise Zr 
is assumed to be CAf{0,a'^), where CAf{0,cF^) denotes the 
circularly symmetric complex Gaussian random variable with 
variance cr^. We assume a block fading scenario, with the 
ratio Hb/Ha denoted as z = je^^, where 7 e M+ and — tt < 
< TT, is referred as the fading state and for simplicity, also 
denoted by (7, 0). Also, it is assumed that z is distributed 
according to a continuous probability distribution. 

Let 5^(^,9) denote the effective constellation at the relay 



during the MA Phase, i.e., 

Snil^O) = {xi + je^^Xj\xi,Xj e S} . 

Let djniniis^^) denote the minimum distance between the 
points in the constellation 5^(7,6'), i.e., 

rfmin{7e^*) = min \{xA-x'^)+'fe^^{xB-x'g)\. (1) 

(^a,=^b),(^'a'^b') 
eSxS 

(^a,^b)¥'(^'a'^'b'> 

From ([T]|, it is clear that there exists values of 76^^ for 
which d„„;„(7eJ') = 0. Let H = {76^^ € C|d™™(7, 6) = 0}. 
The elements of H are said to be the singular fade states. An 
alternative definition of singular fade state is as follows: 

Definition 1: A fade state 76-'^ is said to be a singular fade 
state, if the cardinality of the signal set Suiljd) is less than 
A/2. 

For example, consider the case when symmetric 4-PSK 
signal set used at the nodes A and B, i.e., S = {{±l±j)/y/2}. 
For je^<' = (1 + j)/2, d„unheJ') - 0, since, 

\V2 V2 2 y V2 V2 J ■ 

Alternatively, when 76^^ = {I + j)/2, the constellation 
Sr{j,9) has only 12 (<16) points. Hence 76^" = (1 + j)/2 
is a singular fade state for the case when 4-PSK signal set is 
used at A and B. 

Let {xa,xb) G S X S denote the Maximum Likelihood 
(ML) estimate of {xa, xb) at R based on the received complex 
number Yr, i.e., 

{xa,xb) = arg min |Yr - Hax'a - Hsx'gl (2) 

Broadcast (BC) Phase: Depending on the value of 76^^, R 
chooses a map Ai"^'^ : S x S ^ S' , where S' is the signal set 
(of size between M and AP) used by R during BC phase. 
The elements in 5 x 5 which are mapped on to the same 
complex number in S' by the map M'*-^ are said to form a 
cluster. Let {£1, £2, denote the set of all such clusters. 

The formation of clusters is called clustering, denoted by C^'^. 

The received signals at A and B during the BC phase are 
respectively given by, 

Ya = H'aXr + Za, Yb = H'bXr + Zb, (3) 

where Xr = Ai"''^ {xa, xb) G S' is the complex number 
transmitted by R. The fading coefficients corresponding to the 
R-A and R-B links are denoted by iJ^ and H'g respectively 
and the additive noises Za and Zb are CJV{0, cr^). 

In order to ensure that A (B) is able to decode B's (A's) 
message, the clustering C'^ should satisfy the exclusive law 
0, i.e., 

XT'^(xA,xs) ^ M'''''{x'^,xb), where xa ^ x'^ , e 5, 1 
XT'^(xA,xs) ^ M''-''(xA,x'g), where xb + x'g ,xa&S. J 

(4) 



From an information theoretic perspective, the mapping 
A^T'^ needs to satisfy the exclusive law for the reason outlined 
below. Consider the ideal situation where the additive noises 
at the nodes are zero. It is assumed that the fading state 
76-'^ ^ H. The assumption is required since R can decode 
unambiguously to an element in S x S only if 76-'^ is not 
a singular fade state and is justified since 76-'^ takes values 
from a continuous probability distribution and the cardinality 
of Ti. is finite. During the MA Phase, assume the relay 
jointly decodes correctly to the pair {xa,xb) and transmits 
Xr — M'^'^ {xa,xb) during the BC Phase. The received 
complex symbols at A and B are respectively Ya = H'^^Xr 
and Yb = H'^Xr. At node A, the amount of uncertainty about 
xb which gets resolved after observing Ya, I{xb', Ya/xa) = 
H{xb\xa) - H{xb\Ya,xa) = H{xb) - H{xb\Xr,xa) 
(since xr and xa are independent). Since H{xb\Xr, xa) — 
if and only if the mapping Ai'''^ satisfies the exclusive law, 
the amount of uncertainty about xr (xa) which gets resolved 
at A (B) is maximized if and only if the clustering satisfies 
the exclusive law. 

Definition 2: The cluster distance between a pair of clusters 
£i and Cj is the minimum among all the distances calculated 
between the points xa + jc^^xbjx'^ + ■ye^^x'^ S Sr{j,9), 
where {xa,xb) e A and (x'^jx'^) G Cj. 

Definition 3: The minimum cluster distance of the cluster- 
ing C'^ is the minimum among all the cluster distances, i.e., 

dmin{C~''^) = min | {xa - x'j^) + 76^'* {xb - Xg) |. 

(=^a,xb),(^A'^b) 
esxs, 

MI'O {^A,^B)i^M-''^ 

The minimum cluster distance determines the performance 
during the MA phase of relaying. The performance during 
the BC phase is determined by the minimum distance of the 
signal set S' . Throughout, we restrict ourselves to optimizing 
the performance during the MA phase. For values of 76^^ 
in the neighborhood of the singular fade states, the value 
of dmm{je^^) is greatly reduced, a phenomenon referred as 
distance shortening. To avoid distance shortening, for each 
singular fade state, a clustering needs to be chosen such that 
the minimum cluster distance at the singular fade state is non- 
zero and is also maximized. 

A clustering C^''^ is said to remove a singular fade state h G 
Ti, if the minimum cluster distance of the clustering C^'^^ for 
76-''* = /i is gi-eater than zero. Let Cu = {C^''^ ■.h£'H] de- 
note the set of all such clusterings. Let dmin{C^'^^ 0) denote 
the minimum cluster distance of the clustering C^'''^ evaluated 
at 76-'^. For 76-'^ ^ "H, the clustering C'" is chosen to be C^''^ 

which satisfies d™„(C{''>, 7, 0) > dra^n{C^^'\l,e),yh ^ 

h' G n. 

Note 1: The clusterings which belong to the set C-r need 
not be distinct, since a single clustering can remove more than 
one singular fade state. 

Example 1: In the case of BPSK, if the fade state is 7 = 1 
and 9 = the distance between the pairs (0, 1)(1, 0) is zero as 
in Fig|2];a).The following clustering remove this singular fade 
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(0,1)(1,0) 



BPSK used at node B 
(0,1) 
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(a) 7 = 1 aiid e = (b) 7 = 1 and 9 = jr 



Fig. 2. Effective Constellation at the relay for singular fade states, when the 
end nodes use BPSK constellation. 



State: 

{{(0,1)(1,0)},{(1,1)(0,0)}} 
The minimum cluster distance is non zero for this clustering. 

B. Issues with Koike-Akino Popovski Tarokh's approach 

It is assumed that the channel state information is not 
available at the transmitting nodes A and B during the MA 
phase. A block fading scenario is assumed and the clustering 
used by the relay is indicated to A and B by using overhead 
bits. 

The procedure suggested in ||4J to obtain the set of all 
clusterings, was using a computer search algorithm (closest 
neighbour clustering algorithm), which involved varying the 
fade state values over the entire complex plane, i.e., < 7 < 
00, < 9 < 2tt and finding the clustering for each value 
of channel realization. The total number of network codes 
which would result is known only after the algorithm is run 
for all possible realizations 76^^ which is uncountably infinite 
and hence the number of overhead bits required is not known 
beforehand. Moreover, performing such an exhaustive search is 
extremely difficult in practice, especially when the cardinality 
of the signal set M is large. 

The implementation complexity of the scheme suggested in 
Q is extremely high. It appears that, for each realization of the 
singular fade state, the closest neighbour clustering algorithm 
Q needs to be run at R to find the clustering. In contrast, 
we provide a simple criterion (Section IV A, fTO) ) based on 
which a clustering from the set Cn is chosen depending on 
the value of je^^. In this way, the set of all values of je^^ 
(the complex plane) is quantized into different regions, with a 
clustering from the set C-^ used in a particular region. 

In the closest neighbour clustering algorithm suggested in 
Q, the network coding map is obtained by considering the 
entire distance profile. The disadvantages of such an approach 
are two-fold. 

• Considering the entire distance profile, instead of the 
minimum cluster distance alone which contributes dom- 
inantly to the error probability, results in an extremely 



large number of network coding maps. For example, 
for 8-PSK signal set, the closest neighbour clustering 
algorithm results in more than 5000 maps. 

• The closest neighbour clustering algorithm tries to op- 
timize the entire distance profile, even after clustering 
signal points which contribute the minimum distance. 
As a result, for several channel conditions, the number 
of clusters in the clustering obtained is greater than the 
number of clusters in the clustering obtained by taking the 
minimum distance alone into consideration. This results 
in a degradation in performance during the BC Phase, 
since the relay uses a signal set with cardinality equal to 
the number of clusters. 

In f4\, to overcome the two problems mentioned above, 
another algorithm is proposed, in which for a given 76^^, 
an exhaustive search is performed among all the network 
coding maps obtained using the closest neighbour clustering 
algorithm and a map with minimum number of clusters is 
chosen. The difficulties associated with the implementation of 
the closest neighbour clustering algorithm carry over to the 
implementation of this algorithm as well. 

To avoid all these problems we suggest a scheme, which is 
based on the removal of all the singular fade states. Since the 
number of singular fade states is finite (the exact number of 
singular fade states and their location in the complex plane 
are discussed in Section II), the total number of network 
coding maps used is upper bounded by the number of singular 
fade states. In fact, the total number of network coding maps 
required is shown to be lesser than the total number of singular 
fade states in Section VI. In other words, the total number 
of network coding maps required is known exactly, which 
determines the number of overhead bits required. It is shown 
in section III that the problem of obtaining clusterings which 
remove all the singular fade states reduces to completing a 
finite number of partially filled Latin Squares, which totally 
avoids the problem of performing exhaustive search for an 
uncountably infinite number of values. 

In |4|, an important observation was that for the case of 
QPSK modulation during the MA Phase, there exists several 
channel conditions under which the use of unconventional 
5-ary signal constellation results in a larger throughput. It 
appears as if the mitigation of the multiple access interference, 
comes at the cost of degraded performance during the BC 
phase, because of the use of a signal set with a larger cardi- 
nality during the BC phase. By providing all the clusterings 
explicitly, it is shown in Section III B that the use of a signal 
set with larger cardinality is not required for 8-PSK signal set. 
In other words, the mitigation of multiple access interference 
does not come at the cost of degraded performance during the 
BC phase, for 8-PSK signal set. 

The contributions and organization of the paper are as 
follows: 

• It is shown that the requirement of satisfying the exclusive 
law is same as the clustering being represented by a 
Latin Square and can be used to get the clustering which 
removes singular fade states. In other words, it is shown 



that the problem of finding a clustering which removes 
a singular fade state reduces to filling a partially filled 
Latin Square. 

Using the properties of the set of Latin Squares for a 
given set of parameters, the problem of finding the set of 
maps corresponding to all the singular fade states can be 
simplified to finding the same for only for a small subset 
of singular fade states. Specifically, it is shown that 

1) For the set of all singular fade states lying on a 
circle, from a Latin Square corresponding to one 
singular fade state, Latin Squares for the other 
singular fade states can be obtained by appropriate 
permutation of the columns of the first Latin Square. 

2) There is a one-to-one correspondence between a 
Latin Square corresponding to a singular fade state 
on a circle of radius r and a Latin Square corre- 
sponding to a singular fade state on a circle of radius 
1 

r ' 

It is shown that the bit-wise XOR mapping can remove 
the singular fade state (7 = 1,6' = 0) for any Af-PSK, 
(i.e., for M any power of 2) 

For any Af-PSK signal set, all the clusterings which can 
remove the singularities can be obtained with the aid of 
Latin Squares along with their isotopes. As an example, 
this is shown explicitly for QPSK signal set and 8PSK 
signal set. 

It was shown in |4| that for the case when QPSK signal 
set is used during the MA phase, there exists certain 
values of fade state for which a clustering with cardinality 
5 needs to be used to maximize the minimum cluster 
distance. While the use of clusterings with cardinality 
5, reduces the impact of multiple access interference, 
it adversely impacts the performance during the BC 
phase. For 8-PSK signal set, clusterings which remove 
the singular fade states, all of which have a cardinality 
of 8 are explicitly given, which implies that there is no 
trade-off between the MAC Phase and the BC phase. 
Even though, the problem of completability using M 
symbols of an Af x AI Latin Square in general is 
unsolved, it is shown that the structure of the partially 
filled Latin Square in the problem considered allows the 
construction of explicit Latin Squares, for some singular 
fade states. In other words, by providing some explicit 
constructions, for A/-PSK signal set, it is shown that 
the partially filled Latin Squares corresponding to some 
singular fade states can be completed using M symbols 
itself. 

It was observed in Q that for 4 PSK signal set there 
exists clusterings which remove more than one singular 
fade state. For any Af-PSK signal set, certain special 
cases when multiple singular fade states are removed 
by the same Latin Square are identified. A construction 
algorithm is provided to obtain such Latin Squares. Each 
one of the Latin Squares obtained using the algorithm 
provided removes AP/8 singular fade states. 



• An explicit construction procedure is provided to show 
that certain Latin Squares for Af-PSK signal set, are 
obtainable from the Latin Squares for Af /2-PSK signal 
set. 

• When the end nodes use constellations of different sizes 
Af(= 2^) and iV(= 2''), to get the clusterings it is 
required to fill Latin Rectangles. These Latin rectangles 
are obtained by removing appropriate columns of the 
Latin Squares constructed for the case when both the 
end nodes use constellations of same size which is 
max(Af,7V). 

II. SINGULAR FADE STATES FOR M-PSK SIGNAL 

SET 

Throughout the paper the points in the symmetric Af -PSK 
signal set are assumed to be of the form e-''^'^+^-''^/*^, < 
fc < Af — 1 and Af is of the form 2^, where A is a positive 
integer. Let AS denote the difference constellation of the Af- 
PSK signal set S, i.e., AS = {s, - s'Js,, s'^ € S}. 

For any M-PSK signal set, the set AS is of the form, 

A5 ={0} U |2sin(7rn/M)eJ'=2''/*^|n odd} 

U |2sin(7rn/M)e^('=2'^/"+'^/"V evenj , 

where 1 < n < Af /2 and < fc < Af - L 

In other words, the non-zero points in AS lie on Af/2 
circles of radius 2 sin(7rn/Af ), 1 < n < Af/2 with each circle 
containing Af points. The phase angles of the Af points on 
each circle is 2fc7r/Af, if n is odd and 2fc7r/Af + 7r/Af if n 
is even, where < fc < Af — 1. For example the difference 
constellation for 4-PSK and 8-PSK signal sets are shown in 
Fig. [3] and Fig. [4] respectively. 
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Fig. 3. Difference constellation for 4-PSK signal set 
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Fig. 4. Difference constellation for 8-PSK signal set 



Let us define, 



2 sin(7rn/M)e^'=2'^/*^ if n is odd, 

2 sin(7rn/M)eJ'('=2'^/^+'^/*^) if n is even. 



(5) 



where 1 < n < Af/2 and < fc < M - 1. 

From Q, it follows that the singular fade states are of the 
form. 



7s 



for some Xk,n, Xk' ,n' G 

Lemma 1: For integers fci, fc2, ?i and I2, where 

M 

1 < ki,k2,li,l2 < 7^ and ^ h, 

sin(A;i7r/M) _ sin(ii7r/M) 
sin(fc27r/A'/) ~ smihir/M) ' 

if and only if ki = h and k2 = 

Proof: See fro). ■ 

The following lemma gives the location of the singular fade 
states in the complex plane. 

Lemma 2: The singular fade states other than zero lie on 
M/2+1 circles with M points on each circle, with the 
radii of the circles given by sin(fci7r/M)/ sin(fc27r/M), where 
1 < fci, ^2 < M/2. The phase angles of the M points on each 
one of the circles are given by fc27r/M, < fc < M — 1, if 
both ki and k2 are odd or both are even and kl-K jM -I- 7r/Af, 
< A; < — 1, if only one among k\ and k2 is odd. 

From Lemma [2] it follows that if 76-'^ is a singular fade 
state, ^e^-'^ is also a singular fade state. 

Example 2: For the case when 4-PSK signal set is used 
during the MA Phase, the singular fade states lie on three 
circles as shown in Fig. |5] 



Example 3: For the case when 8-PSK signal set is used 
during the MA Phase, the singular fade states lie on thirteen 
circles as shown in Fig. |6] 
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Fig. 5. Singular fade states for 4-PSK signal set 
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Fig. 6. Singular fade states for 8-PSK signal set 

III. The Exclusive Law and Latin Squares 

Definition 4: Q A Latin Square L of order M on the 
symbols from the set Zj = {0, 1, • • • , i — 1} is an M x M 
array, in which each cell contains one symbol and each symbol 
occurs at most once in each row and column. 

Let the nodes A and B use the same constellation of size Af . 
Consider an M x M array at the relay with the rows (columns) 
indexed by the constellation point used by node A (B), i.e., 
symbols from the set {0, 1, 2, ... , M-\}. The relay is allowed 
to use any constellation with size f > M (using t > M may 
lead to some advantages, see Q). Our aim is to cluster the 
JVP slots in the M x M array such that the exclusive law is 



satisfied. To do so, we will fill in the slots in the array with the 
elements of set Zt in such a way that ( [T5| ) is satisfied, and the 
clusters are obtained by taking all the row-column pairs 
i, j € 0, 1, • • • , Af — 1, such that the entry in the {i, j)— th slot 
is the same symbol from Zf . The specific symbols from Zt are 
not important, but it is the set of clusters that are important. 
Now, it is easy to see that if the exclusive law need to be 
satisfied, then the clustering should be such that an element 
in a row and also in a column cannot be repeated in the same 
row and column. Thus all the relay clusterings which satisfy 
the exclusive law form Latin Squares. Hence, we have the 
following: 

All the relay clusterings which satisfy the exclusive law 
forms Latin Squares, when the end nodes use constellations 
of same size. 

With this observation, the study of clustering which satisfies 
the exclusive law can be equivalently carried out as the study 
of Latin Squares with appropriate parameters. 

A. Removing Singular fade states and Constrained Latin 
Squares 

The relay can manage with constellations of size M in BC 
phase, but it is observed that in some cases relay may not be 
able to remove the singular fade states with t = M and results 
in severe performance degradation in the MA phase |4|. Let 
(fc, l){k\ V) be the pairs which give same point in the effective 
constellation Sr at the relay for a singular fade state, where 
k,k',l,l' e {0,1,. ...,M - 1}. Let fc,fc' be the constellation 
points used by node A and /, be the constellation points used 
by node B. If they are not clustered together, the minimum 
cluster distance will be zero. To avoid this, those pairs should 
be in same cluster This requirement is termed as singularity- 
removal constraint. So, we need to obtain Latin Squares which 
can remove singular fade states and with minimum value for 
t. Therefore, initially we will fill the slots in the M x M 
array such that for the slots corresponding to a singularity- 
removal constraint the same element will be used to fill 
slots. This removes that particular singular fade state. Such 
a partially filled Latin Square is called a Constrained Partial 
Latin Square. After this, to make this a Latin Square, we will 
try to fill the other slots of the partially filled. Constrained 
Partial Latin Square with minimum number of symbols from 
the set Zt. 

Example 4: Consider the case where both A and B uses 
BPSK as shown in Fig. |2] where the effect of noise is 
neglected. There are two singular fade states, one at (7 — 
1,6 = 0) and the other at (7 = 1,6' = tt). We try to 
eliminate the singular fade states one by one. First to remove 
(7 = 1, 6* = 0), the symbol in 0*'' row 1** column (henceforth 
the slot (0,1)) and symbol in (1,0) should be same. Otherwise, 
the minimum cluster distance will be zero. We are using 
symbol 1 (choice of this symbol will not alter the clustering) 
for this and we will get the Constrained Partial Latin Square 
as in Table. [l] This uniquely completes to the Latin Square in 
Table. |ll] Notice that this will remove the singular fade state 
(7 = 1,0 = tt) also. This Latin Square corresponds to the 
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TABLE I 
Partially filled Latin 
Square 



TABLE II 
Completely Filled Latin 
Square 



bit-wise XOR mapping, but with higher order constellations 
the number of singular fade states increases and bit-wise XOR 
cannot remove (will be seen in the sequel) all the singular fade 
states. 

Let the M-PSK points be e^i2k+i)7r/M ^ ^^^^.^ < k < 
M — 1. For simplicity, by the point k we mean the point 
gj(2/c+i)7r/M rpj^g following lemma shows that bit-wise XOR 
mapping removes the singular fade state (7 = 1,9 = 0), for 
any i\f-PSK signal set. 

Lemma 3: When the user nodes use 2'^-PSK constellations, 
the singular fade state (7 = 1,0 = 0) is removed by bit-wise 
XOR mapping (denoted by ©), for all A. 

Proof: For {{k,k'), {1,1')} to be a singularity constraint 
corresponding to the singular fade state (7 = 1,6 — 0), 



e M 



e M 



■ e M 



sin [7r(fc - k')/2^] 
sin [7r(F - /)/2^] 



1, i.e.. 



J M = 1. 



(6) 



Equating the amplitudes on both sides of (|6]l, one gets 

sin [nik - k')/2^] = sin [Tr{l' - l)/2^] , 



which leads to the following two cases: 
Case ( i): k — k' = I' — I 



Cases ( ii): 



TT{k-k') 



nil' -I) 



k~l = k'-l'+2^ 



2^ 2^ _, 

Equating the phases on both the sides of Q leads to 



2^ 



{k + k' - I - I') ^ Q =^ k + k' ^l + l'. 



(7) 



Combining Case (i) above and (|7]i gives {k',l') = {l,k), i.e., 
the singularity-removal constraint is of the form {(fc, l){l, k)}. 
In other words, the clustering should satisfy this symmetry. 

Combining Case (ii) above and (|7]i leads to k = I + 2^~^ 
irrespective of k' , I' and k' — I' + 2^~^ irrespective of k, I. In 
other words, {{I + 2^-^,1)}, I e {0, 1, 2^} is the set of 
singularity -removal constraints. 

From the above, one can conclude that a clustering which 
removes the singular fade state (7 = 1,0 = 0) should have 

(i) A symmetric Latin Square, meaning that the cells (fc, I) 
and {I, k) should have the same symbol. 

(ii) A Latin Square with the symbols in the cells {{I + 
2^"\ /)}, and I e {0, 1, ....,2"^} being the same. 

The Latin Square produced by bit-wise XOR mapping is 
clearly symmetric. Moreover, the quantity {I + 2^^^) / is 
always equal to 2^^^ for all values of I, i.e., the symbols in 
all the cells of the set {{1 + 2^-^,1)}, I e {0, 1, 2^} are 
the same. Hence the XOR map removes the singular fade state 
(7 = 1,0 = 0). ■ 



Definition 5: |8| Two Latin Squares L and L ' (using the same 
symbol set) are isotopic if there is a triple (f,g,h), where / is a 
row permutation, 5 is a column permutation and his a symbol 
permutation, such that applying these permutations on L gives 
L'. 

Lemma 4: Two Latin Squares L and L' which remove the 
singular fade states (7, 0) and (7,6*'), respectively, (i.e., two 
singular fade states on the same circle), are Isotopic that are 
obtainable one from another by a column permutation alone. 
If 9' — 9 = kjj, L' can be obtained by cyclic shifting of the 
columns of L, k times in the anticlockwise direction. 

Proof: Let L and L', respectively remove the singular 
fade states (7,6*) and (7,6*'). 

The effect of rotation in the plane by an angle 9' — 9 
due to channel fade coefficients Ha and Hb can be viewed 
equivalently as a relative rotation of the constellation used 
by B by an angle 9' — 9 with respect to the constellation 
used by A and no relative rotation between the channel fade 
coefficients and Hb ■ Let S and S' be the resulting rotated 
constellations after rotation in the constellation of node B 
corresponding to an angle 9' — 9. 

Since there are M singular fade states for a specific 7, 



(shown in |10|), and they are all spaced by same angular 
separation, 9' — 9 is an integer multiple of 2tt/M which is an 
angular separation of the i\/-PSK constellation points. That is, 
a rotation in the channel by an angle 6*' — 6' is equivalent to a 
rotation in the constellation points in the M -PSK constellation. 
So, we can obtain the Latin Square L' by column permutations 
in L, since the columns are indexed by constellation points 
used by node B. This means, if we obtain the Latin Square 
for a singular fade state (7,6*), then by appropriately shifting 
the columns we obtain the Latin Squares that remove all the 
other singular fade states of the form (7,0')- This completes 
the proof. ■ 

From |4] it follows that for each circle, it is enough if we 
obtain one Latin Square which removes a singular fade state 
on that circle. The Latin Squares which remove the other 
singular fade states can be obtained by column permutation. 
For example, all the Latin Squares which remove the singular 
fade states on the unit circle can be obtained from the bit-wise 
XOR map by column permutation. 

Definition 6: A Latin Square is said to be the Transpose 
of a Latin Square L, if L^[i,j) = L{j,i) for all i,j G 
{0,1,2,..,A/-1}. 

Recall from Section II that if ^e^^ is a singular fade state. 



then is also a singular fade state. The following Lemma 



shows that the transpose of the Latin Square which removes 
76^^ removes the singular fade state ■^e~^^ . 

Lemma 5: If the Latin Square L removes the singular fade 
state (7,6*), then the Latin Square removes the singular 
fade state {^,-9). 

Proof: Let {(fci, /i)(fc2, ^2)} be a singularity-removal 
constraint for the singular fade state (7, 0). Then, 



1 


1 





2 


3 


2 


3 



QPSK used at node A 



QPSK used at node B 



Fig. 7. QPSK Constellations used at the end nodes 



sm 



/ 7r(fci-fc2) \ 
V M ) 



M 



and 9 



M 



(fci 



Taking transpose in the constraint we will obtain 
{(^1, fci)(^2, ^2}. Let this constraint correspond to the 
singular fade state {j',9'). Then, 

^ sin(^) sin(^) 

T = — = — = lH 

sin(^%=M) sin(^^(^) 

Similarly, 

9' = ^{l, + l^-ki- fca) = -^Aki + k2-li- h) = -e. 
Ivl Ivl 

This completes the proof. ■ 
Lemma |5] implies that it is enough to obtain the Latin 
Squares which remove those singular fade states which lie 
on or inside the unit circle centered at the origin. The Latin 
Squares which remove the singular fade states which lie 
outside the unit circle can be obtained by taking transpose. 

IV. Illustrations 

A. End nodes use QPSK 

There are 12 singular fade states (shown in Fig|5]l, when 
both the end nodes A and B use QPSK as in Fig|7] 
The singular fade states are 

7 = 1; e' = 0,+7r/2,-7r/2,7r 
7=1/^/2; = +7r/4,+37r/4, -7r/4,-37r/4 
7 = ^2; 6* = +7r/4,+37r/4, -7r/4,-37r/4. 

We remove singular fade states one by one. Consider first, 
the case (7 = 1,0 = 0). The singularity -removal constraints 
are 

{(0,1)(1,0)}; {(0,2)(1,3)(2,0)(3,1)}; {(0, 3)(3, 0)}; 
{(1,2)(2,1)}; {(2,3)(3,2)}. 

Satisfying these constraints, a Latin Square can be con- 
structed with i=4, in three different ways, Li,L2 and L3 



as shown in figure 10(a) figure 10(b) and figure 10(c) All 
these three clusterings corresponding to each Latin Square give 
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Ll With symbols permuted. 

Fig. 8. Two rotations of Li gives same clustering as that of Li 



Fig. 9. Obtaining Le from L\ by column shifting 



the same performance on the basis of first minimum cluster 
distance in the MA phase. But the advantage with the one 



shown in figure 10(a) is that it removes singular fade state at 
(7 — 1,6* = tt) also. This is easily verified, by seeing that 
after two cyclic shifts in the columns of Li the clustering that 
it results in is the same as the old one. This is explicitly shown 
inFig|8] 

The singularity-removal constraints for (7 = 1,6 — it) are 

{(0,3)(1,2)}; {(0,0)(1,1)(2,2)(3,3)}; {(0, 1)(3, 2)}; 
{(1,0)(2,3)}; {(2,1)(3,0)}. 

The Latin squares to remove this singular fade state, ^1,^4 



and L5 are shown in figure 10(a) figure 11(a) and figure 11(b) 
respectively. 

In order to reduce the total number of different clusterings 
we select the clustering corresponds to the Latin Square Li 



shown in figure 10(a) as the clustering to remove both these 
singular fade states. The corresponding clustering, Co is 



{(0,1)(1,0)(2,3)(3,2)}, 
{(0,3)(3,0)(1,2)(2,1)}, 



{(0,2)(1,3)(2,0)(3,1)}, 
{(0,0)(1,1)(2,2)(3,3)}. 



Now consider the singular fade state (7 = 1,0 ~ '^l'^)- The 
singularity-removal constraints are 

{(0,0)(1,3)}; {(0,1)(1,2)(2,3)(3,0)}; {(0,2)(3,3)} 
{(1,1)(2,0)}; {(2,2)(3,1)} 

In this case also Latin Square can be constructed with t ~ 4, 
in three different ways as shown in figure |12(a)| figure |12(b)| 
and figure 12(c) But as in earlier case out of these three one. 



Lg (shown in figure |12(a) 1 will remove singular fade state 
(7 = 1,9 = — 7r/2). The singularity-removal constraints for 

(7 = 1,61 = -7r/2) are 

{(0,2)(1,1)}; {(0,3)(1,0)(2,1)(3,2)}; {(0, 0)(3, 1)}; 
{(1,3)(2,2)}; {(2,0)(3,3)}. 

All the Latin Squares which remove the singular fade state 
(7 = 1,6* = —tt/2) are shown in figure 12(a) figure 13(a) and 
figure 13(b) We will select that clustering which reduces total 
number of different clusterings, i.e, Lq as was done before. 
The corresponding clustering, Ci is as follows: 



{(0,0)(1,3)(2,2)(3,1)}, 
{(0,2)(3,3)(1,1)(2,0)}, 



{(0,1)(1,2)(2,3)(3,0)}, 
{(0,3)(1,0)(2,1)(3,2)}. 



Lg are Isotopic Latin Squares. That is, clustering correspond- 
ing to Lq is obtained by, cyclically shifting the columns of 
Ll (since columns are indexed by constellation points used 
by node B) as in Fig|9] 

Now we have removed four singular fade states till now all 
of them on the unit circle. Consider next, (7 = 1/^/2,9 = 
7r/4). The singularity -removal constraints are 

{(0,1)(1,3)}; {(0,2)(3,0)}; {(1, 2)(2, 0)}; {(2, 3)(3, 1)}. 

The corresponding partially filled Latin Square is shown in 



figure 14(a) 



It cannot be completed with t=4. That means that the relay 
has to use a constellation of size more than four. We can see 
that it can be completed with t=5. We get two clusterings as 
given in figure |15(a)| and figure |15(b)| Both will remove this 
singular fade state and use constellation of size five. 

The clustering corresponding to Latin Square Lu, denoted 
by C2 is 

{(0,0)(2,3)(3,1)}, {(0,1)(1,3)(2,2)}, {(0,2)(1,1)(3,0)}, 
{(0,3)(1,0)(2,1)(3,2)}, {(1,2)(2,0)(3,3)}. 

The clustering corresponding to Latin Square L12, denoted 
by C3 is 

{(0,0)(1,1)(2,2)(3,3)}, {(0,1)(1,3)(3,2)}, {(0,2)(2,1)(3,0)}, 
{(0,3)(1,2)(2,0)}, {(1,0)(2,3)(3,1)}}. 

Now considering the next singular fade state (7 = y/2, 9 = 
7r/4), by the same procedure as before, the singularity-removal 
constraints are 

{(0,1)(2,0)}, {(0,2)(2,3)}, {(1,2)(3,1)}, {(1,3)(3,0)}. 

The partially filled Latin Square is shown in figure [T4(b)l This 
cannot be completed with t = 4, but by i = 5 it can be 
completed in two ways as in L13 and L14 shown in figure 16(a) 



and figure 16(b) The corresponding clusterings are shown in 
the Table [ml 
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Fig. 10. Completely Filled Latin Squares for 7 = 1,6 = 



The interesting point here is that the Latin Squares Li and 



TABLE III 

Clusterings Obtained for different singular fade states when the end nodes use QPSK constellations 



SI No 


Singular fsdc st3.t6s 


Clustering 


Cluster 


1 


7 = 1,61 = 


Co 


{{(0,1)(1,0)(2,3)(3,2)},{(0,2)(1,3)(2,0)(3,1)},{(0,3)(3,0)(1,2)(2,1)},{(0,0)(1,1)(2,2)(3,3)}} 


2 


— 1 6 — -K 11 

1 — -L, 1^ — n 1 ^ 




f-fffl 0¥1 3V2 2in 1^1^ UO lin 21^2 3V3 011^ UO 2V3 3V1 11C2 0^1^ UO 3¥1 OIQ 1¥3 2111 


3 


-y = 1,6» = TT 


Co 


{{(0,1)(1,0)(2,3)(3,2)},{(0,2)(1,3)(2,0)(3,1)},{(0,3)(3,0)(1,2)(2,1)},{(0,0)(1,1)(2,2)(3,3)}} 


4 


7 = 1, e = -7r/2 


Ci 


11(0 0)(1 3K2 2)(3 1)>,I(0 1)(1 2)(2 3)(3 0)>,l(0 2)(3 3)(1 1)(2 0)>,l(0 3)(1 0)(2 1)(3 2)11 


5.a 
5.b 


7 = 1/^2,6* = 7r/4 


C2 

C3 


{{(0,0)(2,3)(3,1)},{(0,1)(1,3)(2,2)},{(0,2)(1,1)(3,0)},{(0,3)(1,0)(2,1)(3,2)},{(1,2)(2,0)(3,3)}} 
||(0,0)(1,1)(2,2)(3,3)>,|(0,1)(1,3)(3,2)1,|(0,2)(2,1)(3,0)>,|(0,3)(1,2)(2,0)>, 1(1,0) 


6.a 
6.b 


7 = v^,6» = 7r/4 


Ci 
C5 


{{(0,0)(1,1)(2,2)(3,3)},{(0,1)(2,0)(3,2)},{(0,2)(1,0)(2,3)},{(0,3)(1,2)(3,1)},{(1,3)(2,1)(3,0)}} 
11(0 0)(1 2)(3 1)>,I(0 1)(2 0)(3 3)>,l(0 2)(1 1)(2 3)>,l(0 3)(1 0)(2 1)(3 2)1,1(1 3)(2 2)(3 0)U 


7.a 
7.b 


7 = l/v^,6l = 37r/4 


Cs 
C5 


{{(0,0)(1,2)(2,1)},{(0,1)(1,0)(3,3)},{(0,2)(1,3)(2,0)(3,1)},{(0,3)(2,2)(3,0)},{(1,1)(2,3)(3,2)}} 
11(0 0)(1 2)(3 1)},((0 1)(2 0)(3 3)>,l(0 2)(1 1)(2 3)>,l(0 3)(1 0)(2 1)(3 2)1,1(1 3)(2 2)(3 0)>> 


8.a 
8.b 


7 = v^,e = 37r/4 


C2 
C7 


{{(0,0)(2,3)(3,1)},{(0,1)(1,3)(2,2)},{(0,2)(1,1)(3,0)},{(0,3)(1,0)(2,1)(3,2)},{(1,2)(2,0)(3,3)}} 
{{(0,0)(2,3)(3,2)},{(0,1)(1,0)(2,2)},{(0,2)(1,3)(2,0)(3,1)},{(0,3)(1,1)(3,0)},{(1,2)(2,1)(3,3)}} 


9.a 
9.b 


7 = 1/^2, e = -37r/4 


C7 
Cg 


{{(0,0)(2,3)(3,2)},{(0,1)(1,0)(2,2)},{(0,2)(1,3)(2,0)(3,1)},{(0,3)(1,1)(3,0)},{(1,2)(2,1)(3,3)}} 
{{(0,0)(1,3)(3,2)},{(0,1)(1,2)(2,3)(3,0)},{(0,2)(2,1)(3,3)},{(0,3)(1,1)(2,0)},{(1,0)(2,2)(3,1)}} 


lO.a 
lO.b 


7 = v^,6» = -37r/4 


Cs 
C9 


{{(0,0)(1,2)(2,1)},{(0,1)(1,0)(3,3)},{(0,2)(1,3)(2,0)(3,1)},{(0,3)(2,2)(3,0)},{(1,1)(2,3)(3,2)}} 
{{(0,0)(1,3)(2,1)},{(0,1)(1,2)(2,3)(3,0)},{(0,2)(1,0)(3,3)},{(0,3)(2,2)(3,1)},{(1,1)(2,0)(3,2)}} 


U.a 
ll.b 


7 = l/v^,0 = -7r/4 


C4 
C9 


{{(0,0)(1,1)(2,2)(3,3)},{(0,1)(2,0)(3,2)},{(0,2)(1,0)(2,3)},{(0,3)(1,2)(3,1)},{(1,3)(2,1)(3,0)}} 
{{(0,0)(1,3)(2,I)},{(0,1)(1,2)(2,3)(3,0)},{(0,2)(1,0)(3,3)},{(0,3)(2,2)(3,1)},{(1,1)(2,0)(3,2)}} 


12.a 
12.b 


7 = v^,e = -7r/4 


C3 
Cg 


{{(0,0)(1,1)(2,2)(3,3)},{(0,1)(1,3)(3,2)},{(0,2)(2,1)(3,0)},{(0,3)(1,2)(2,0)},{(1,0)(2,3)(3,1)}} 
{{(0,0)(1,3)(3,2)},{(0,1)(1,2)(2,3)(3,0)},{(0,2)(2,1)(3,3)},{(0,3)(1,1)(2,0)},{(1,0)(2,2)(3,1)}} 
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Fig. 11. Completely Filled Latin Squares for 7 = 1, 



Fig. 13. Completely Filled Latin Squares for 7 = 1, 
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Fig. 1 2. Completely Filled Latin Squares for -y = 1,6 = tt/2 

Next, consider the singular fade state {-y — 1/ \/2, 9 ~ 
37r/4). The singularity -removal constraints are 

{(0,0)(1,2)}, {(0,1)(3,3)}, {(1,1)(2,3)}, {(2,2)(3,0)}. 

The corresponding Latin Squares are shown in figure |16(b)| 

V2,0 



and figure 17(a) 
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Fig. 14. Partially Filled Latin Squares for 7 
7 = ^,0 = 7r/4 



I/V2, 8 = 7r/4 and for 



(7 = 72,6* = -37r/4) are 

{(0,0)(2,1)}, {(0,3X2,2)}, {(1,0)(3,3)}, {(1,1)(3,2)}. 



and the Latin Squares are given in figure |17(a)| and figure 
18(b) The singularity -removal constraints for singular fade 



For the singular fade state (7 
singularity-removal constraints are 



37r/4), the state (7 = l/\/2, B = -7r/4) are 

{(0,2)(1,0)}, {(0,3)(3,1)}, {(1,3)(2,1)}, {(2,0)(3,2)}. 



{(0,0)(2,3)}, {(0,1)(2,2)}, {(1,1)(3,0)}, {(1,2)(3,3)} 



and the corresponding Latin Squares are shown in figure [T5(a)| 
and figure |17(b)| . 

Similarly, for the singular fade state (7 = 1/ a/2, = 
— 37r/4), the singularity-removal constraints are 

{(0,0)(3,2)}, {(0,3)(1,1)}, {(1,0)(2,2)}, {(2,1)(3,3)} 



The Latin Squares are given in figure fT6(a)| and figure [18(5)] 
The singularity-removal constraints for singular fade state 7 = 
\/2 and 6 = -ir/A are 

{(0,2)(2,1)}, {(0,3)(2,0)}, {(1,0)(3,1)}, {(1,3)(3,2)}. 



The Latin Squares are given in figure 15(b)| and figure |18(a) 



with the corresponding Latin Squares as shown in figure [T7(b) 
and figure |18(a)| 

The singularity-removal constraints for singular fade state 



It is observed that to remove all other singular fade states not 
lying on unit circle the relay needs a constellation of size five. 



Table III shows the singular fade states and the corresponding 
clusterings. There are two clusterings to remove a singular fade 
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Fig. 15. Lii for 7 = 1/^2, 9 = tt/A and 7 = ^2, 6» 
for 7 = 1 /v^, e = 7r/4 and 7 = \/2, 6> = -7r/4 






1 


2 


3 


1 


4 





2 


2 





3 


4 


3 


2 


4 


1 






1 


2 


3 


1 


3 


4 


2 


2 


4 


1 





3 


2 





4 



= 37r/4 and L12 Fig. 17. L15 for 7 = v^, 6 = -37r/4 and 7 = 1/^2, 6» = 37r/4 and Lie 
V2, 9 = 37r/4 and 7 = l/v/2, 6» = -37r/4 



for 7 : 



4 





1 


2 


1 


4 


2 


3 





3 


4 


1 


3 


2 





4 



2 





1 


4 


4 


1 


2 


3 





4 


3 


1 


3 


2 


4 






"(a 

Fig. 16. Li3 for 7 = V2, 6» = 7r/4 and 7 = 1/^2, 9 = -7r/4 and L14 
' = ^2,9 = 7r/4 and 7 = 1/a/2,6I = 37r/4 



for 7 : 
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Fig. 18. Li7 for 7 = v^, 6» = -7r/4 and 7 = 1/v^, 6 = -37r/4 and 
Lis for 7 = V2, 6» = -37r/4 and 7 = 1/^2, d = -7r/4 



State for all singular fade states except for those with 7 = 1. 
We can select any one. Anyone from the two {C2:Cz] can be 
selected to remove singular fade state (7 = l/y/2,9 = 7r/4). 
After that, by column permutations we can remove the singular 
fade states with (7 = 1/^2, and 9 = +37r/4, -7r/4, -37r/4. 
By taking transpose of the Latin Square for (7 = l/\/2, 9 = 
7r/4) we can remove singular fade state (7 = \/2, 9 = —i: /A). 
After that by column permutations we can remove the singular 
fade states with 7 = \/2 and 61 = +37r/4, +7r/4, -37r/4. If 
we select C2 to remove (7 ~ l/\/2, 6* = 7r/4), we will get the 
following set of clusterings {Co;Ci,C2,C4,C6,C8} to remove 
all the singular fade states. In the other case, when we select C3 
to remove (7 = 1/a/2, 6* = 7r/4) we will get the following set 
of clusterings {Cq , Ci , C3 , C5 , C7 , Cg } to remove all the singular 
fade states. 

B. End nodes use 8PSK 

Consider the scenario where the end nodes use 8-PSK 



sin7r/8 sin7r/8 sin27r/8 sin37r/8 sin47r/8 



constellation as shown in Fig 19 



For 8-PSK signal set, the 104 singular fade states are 
shown in Fig|6] These singular fade states lie in 13 circles, 
8 in each circle separated by an angle tt/4. It is sufficient 
to get Latin Squares for these 13 values, then by isotopic 
property of the Latin Squares all other singular fade states are 
removed. Further by taking transposes of the Latin Square for 
a singular fade state 7,6*, another singular fade state I/7, —9 
is removed. As discussed XOR removes singular fade state 
7 = 1,6' = 0. Finally it is seen that obtaining six Latin 
Squares for different singular fade state values for 7 < 1 (or 
7 > 1) and the XOR is enough to remove all the 104 singular 
fade states. The 104 singular fade states are, 
for 9 — 2m7r/8 , where m = 0, 1, 2 • • • 7 



smTT/ 



sm TT / 



sin37r/ 



7 



sin7r/8' sin37r/8' sin7r/8 ' sin27r/8 
for 9 = {2m + l)7r/8 , where m = 0, 1, 2 • 



sin47r/8 sin27r/ 
and 



sin 47r/8' 



sin47r/8' sin27r/8' sin37r/8' sin47r/8' sin37r/8' 

sin37r/8 sin27r/8 , sin47r/8 

— , — and — . 

sin27r/8 sinvr/S sin7r/8 
We consider singular fade states one by one. The singularity 

constraints for 7 = 1,9 = are 
{(0,1)(1,0)}, {(0,2)(2,0)}, {(0,3)(3,0)}, {(0,5)(5,0)} 



{(0,6)(6,0)}, 
{(1,4)(4,1)}, 
{(2,4)(4,2)}, 
{(3,5)(5,3)}, 
{(4,7)(7,4)}, 



{(0,7)(7,0)}, 
{(1,6)(6,1)}, 
{(2,6)(6,2)}, 
{(3,6)(6,3)}, 
{(5,6)(6,5)}, 



{(1,2)(2,1)}, 
{(1,7)(7,1)}, 
{(2,7)(7,2)}, 
{(4,5)(5,4)}, 
{(5,7)(7,5)}, 



{(1,3)(3,1)} 
{(2,3)(3,2)} 
{(3,4)(4,3)} 
{(4,6)(6,4)} 
{(6,7)(7,6)} 



{(0, 4)(1, 5)(2, 6)(3, 7)(4, 0)(5, 1)(6, 2)(7, 3)} 

XOR can remove this singular fade state. Further by ap- 
propriately column shifting all the singular fade states with 
7 = 1 is removed. For example consider 7 = 1, = 7r/4. The 
singularity constraints are 



{(0,0)(1,7)}, 
{(0,5)(6,7)}, 
{(1,3)(4,0)}, 
{(2,3)(4,1)}, 



{(0,1)(2,7)}, 
{(0,6)(7,7)}, 
{(1,5)(6,0)}, 
{(2,5)(6,1)}, 



{(0,2)(3,7)}, 
{(1,1)(2,0)}, 
{(1,6)(7,0)}, 
{(2,6)(7,1)}, 



{(0,4)(5,7)} 
{(1,2)(3,0)} 
{(2,2)(3,1)} 
{(3,3)(4,2)} 



{(3,4)(5,2)}, {(3,5)(6,2)}, {(4, 4)(5, 3)}, {(4,5)(6,3)} 
{(4,6)(7,3)}, {(5,5)(6,4)}, {(5,6)(7,4)}, {(6,6)(7,5)} 
{(0, 3)(1,4)(2, 5)(3, 6)(4, 7)(5, 0)(6, 1)(7, 2)} 




Fig. 19. 8PSK constellation used by nodes A and B 



Next consider 7 
singularity constraints are 



sin Stt/S 



0.414,61 



0. The 



{(0,2)(1,7)}, {(0,3)(1,6)}, {(0,5)(7,2)}, {(0,6)(7,1)} 
{(1,3)(2,0)}, {(1,4)(2,7)}, {(2,4)(3,1)}, {(2,5)(3,0)} 
{(3,5)(4,2)}, {(3,6)(4,1)}, {(4, 6)(5, 3)}, {(4,7)(5,2)} 
{(5,0)(6,3)}, {(5,7)(6,4)}, {(6, 0)(7, 5)}, {(6,1)(7,4)} 
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sinvr/S 


fl — 









sin Svr/S 
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(d) 7: 



-,e = 
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(a) 7 



sinvr/S ' 



7r/4 



sin27r/8 
(b) 7 = — r^,0 = O 
sm47r/8 



sin27r/8 
(c) 7 = — ^,e = 7r/2 
sin 471/8 
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sin 37r/8 ' 

Fig. 20. Latin Squares Corresponding to Different singular fade states 



The Latin Square shown in figure 20(c) removes this singu- 
lar fade state. Further by column shifting all the singular fade 

states with 7 = 0.414 are removed. Consider 7 = 

0.414, 6* = 7r/4. The singularity constraints are 



sin 37r/8 



sin47r/8 
(d) 7= . ^ ' ,e = -K/2 
sin 27r/8 

Fig. 21. Latin Squares CoiTesponding to Different singular fade states 



{(0,1)(1,6)}, {(0,2)(1,5)}, {(0,4)(7,1)}, {(0,5)(7,0)} 
{(1,2)(2,7)}, {(1,3)(2,6)}, {(2,3)(3,0)}, {(2,4)(3,7)} 
{(3,4)(4,1)}, {(3,5)(4,0)}, {(4,5)(5,2)}, {(4,6)(5,1)} 
{(5,7)(6,2)}, {(5,6)(6,3)}, {(6, 7)(7, 4)}, {(6,0)(7,3)} 

The Latin Square which removes the above singular fade state 
is shown in figure |20(d)] Latin Squares which removes the 

singular fade states with 7 — — — are obtained by taking 

sm7r/8 

transposes of the Latin Squares which removes the singular 



TABLE IV 

Clusterings Obtained for different singular fade states on each circle (for 7 < 1) when the end nodes use 8-PSK constellations 



Sl.No 



Singular fade states 



Cluster 



{{(0, 1)(1, 0)(2, 3)(3, 2)(5, 4){4, 5)(6, 7)(7, 6)}, {(0, 2)(1, 3)(2, 0)(3, 1)(4, 6)(6, 4)(5, 7)(7, 5)} 
{{0, 3)(3, 0)(1, 2)(2, 1)(4, 7)(5, 6)(6, 5){7, 4)}, {(0, 4)(4, 0)(1, 5)(5, 1)(3, 7)(2, 6)(6, 2)(7, 3)} 
{{0, 5)(1, 4){2, 7)(3, 6)(4, 1)(5, 0)(6, 3){7, 2)}, {(0, 6)(6, 0)(1, 7)(7, 1)(4, 2)(5, 3)(3, 5)(2, 4)} 

{(0, 7)(7, 0)(1, 6)(6, 1)(4, 3)(5, 2)(2, 5)(3, 4)}, {(0, 0)(1, 1)(2, 2)(3, 3)(4, 4)(5, 5)(6, 6), (7, 7)}} 



7 = i,e = 



{{(0, 0)(1, 5)(2, 2)(3, 7)(4, 4)(5, 1)(6, 6)(7, 3)}, |(0, 1)(1, 4){2, 7)(3, 2)(4, 5)(5, 0)(6, 6){7, 3)} 
{(0, 2)(3, 7)(2, 4)(4, 6)(5, 3)(6, 0){7, 5)}, {(0, 3)(1, 6)(2, 1)(3, 4){4, 7)(5, 2)(6, 5)(7, 0)} 

{(0, 4)(1, 1)(2, 6)(3, 3)(4, 0)(5, 5)(6, 2)(7, 7)}, {(0, 5)(3, 6)(1, 0)(2, 3)(4, 1)(5, 4)(6, 7)(7, 2)} 

{{0, 6)(3, 5)(1, 3)(2, 0)(4, 2){5, 3)(6, 4)(7, 1)},{(0, 7)(3, 0)(1, 2){2, 5)(4, 3){5, 6)(6, 1)(7, 4)}} 



sm 7r/8 

7 = — — = 



sin 37r/8 



11(0, 0)(1, 7)(2, 1)(3, 2)(4, 5)(5, 6)(6, 4)(7, 3)|, |(0, 1)(1, 2)(2, 0)(3, 3)(4, 4)(5, 7)(6, 5)(7, 6)| 
{(0, 2)(3, 1)(1, 5)(2, 6)(4, 3)(5, 4)(6, 7)(7, 0)}, {(0, 3)(1, 4)(2, 7)(3, 0)(4, 6)(5, 5)(6, 2)(7, 1)} 

{(0, 4)(1, 3){2, 5)(3, 6)(4, 1)(5, 2)(6, 0)(7, 7)}, {(0, 5)(1, 6)(3, 7)(2, 4){4, 0)(5, 3)(6, 1)(7, 2)} 

{(0, 6)(3, 5)(1, 1)(2, 2)(4, 7)(5, 0)(6, 3)(7, 4)}.{(0, 7)(1, 0)(2, 3)(3, 4)(4, 2)(5, 1)(6, 6)(7, 5)}} 
{{{0, 0)(1, 4)(2, 6)(3, 5)(4, 1){5, 7)(6, 3)(7, 2)}, {(0, 1)(1, 2)(2, 4)(3, 0)(4, 5)(5, 6)(6, 7){7, 3)} 
{(0, 2)(1, 6){2, 5)(3, 4)(4, 0)(5, 3)(6, 1)(7, 7)}, {(0, 3)(1, 7)(2, 1)(3, 2)(4, 4)(5, 0)(6, 5)(7, 6)} 

{{0, 4)(1, 3){2, 7)(3, 1)(4, 2)(5, 5)(6, 6)(7, 0)}, {(0, 5)(1, 0)(2, 3)(3, 6)(4, 7)(5, 4)(6, 2)(7, 1)} 

{(0, 6)(1, 5){2, 2)(3, 7){4, 3)(5, 1)(6, 0)(7, 4)}, {(0, 7){1, 1)(2, 0)(3, 3)(4, 6)(5, 2)(6, 4)(7, 5)}} 



sin 27r/8 

7 = -r^& = 

sin47r/8 



sin 7r/8 
sin47r/8' 



{{(0, 0)(1, 6)(2, 5)(3, 7)(4, 4){5, 2)(6, 1)(7, 3)}, {(0, 1)(1, 7)(2, 2)(3, 0)(4, 5)(5, 3)(6, 6)(7, 4)} 
{(0, 2)(1, 0){2, 3)(3, 1)(4, 6)(5, 4)(6, 7)(7, 5)}, {(0, 3)(1, 5)(2, 4)(3, 2){4, 7)(5, 1)(6, 0)(7, 6)} 

{(0, 4)(1, 3){2, 1)(3, 6)(4, 0)(5, 7)(6, 5){7, 2)}, {(0, 5)(1, 2)(2, 0)(3, 3){4, 1)(5, 6)(6, 4)(7, 7)} 

{(0, 6)(1, 1)(2, 7)(3, 4)(4, 2)(5, 5)(6, 3)(7, 0)}, {(0, 7)(1, 4)(2, 6)(3, 5)(4, 3)(5, 0)(6, 2)(7, 1)}} 
{{{0, 0)(1, 6)(2, 3)(3, 5)(4, 4){5, 2)(6, 7)(7, 1)}, {(0, 1)(1, 4)(2, 2)(3, 7)(4, 5)(5, 0)(6, 6){7, 3)} 
{{0, 2)(1, 0)(2, 7)(3, 1)(4, 6)(5, 4)(6, 3){7, 5)}, {(0, 3)(1, 1)(2, 6)(3, 4)(4, 7)(5, 5)(6, 2)(7, 0)} 

{(0, 4)(1, 7)(2, 5)(3, 2)(4, 0)(5, 3)(6, 1)(7, 6)}, {(0, 5)(1, 2)(2, 4)(3, 3)(4, 1)(5, 6)(6, 0)(7, 7)} 

{(0, 6)(1, 3){2, 1)(3, 0)(4, 2)(5, 7)(6, 5)(7, 4)}, {(0, 7){1, 5)(2, 0)(3, 6)(4, 3)(5, 1)(6, 4)(7, 2)}} 



sin 27r/8 



sin 27r/8 
sin 37r/8 



{{(0, 0)(1, 2)(2, 3)(3, 4)(4, 1)(5, 7)(6, 5)(7, 6)}, {(0, 1)(1, 7)(2, 6)(3, 0)(4, 2)(5, 5)(6, 3)(7, 4)} 
{(0, 2)(1, 6){2, 7)(3, 5)(4, 4)(5, 3)(6, 0){7, 1)}, {(0, 3)(1, 4)(2, 5)(3, 2){4, 0)(5, 1)(6, 6)(7, 7)} 

{(0, 4)(1, 0)(2, 1)(3, 3)(4, 6)(5, 2)(6, 7){7, 5)}, {(0, 5)(1, 1)(2, 0)(3, 6){4, 7)(5, 4)(6, 2)(7, 3)} 

{(0, 6)(1, 5)(2, 4)(3, 7)(4, 3)(5, 0)(6, 7)(7, 2)}, {(0, 7)(1, 3)(2, 2)(3, 1)(4, 5)(5, 6)(6, 4)(7, 0)}} 



sin37r/8 ^ 
sin47r/8' 



7r/8 



fade states with 7 
sin37r/8 



For example, consider 7 



sinSTT/ 

7r/4. The singularity constraints are 



{(1,0)(6,1)}, {(2,0)(5,1)}, {(4,0)(1,7)}, {(5,0)(0,7)} 
{(2,1)(7,2)}, {(3,1)(6,2)}, {(3,2)(0,3)}, {(4,2)(7,3)} 
{(4,3)(1,4)}, {(5,3)(0,4)}, {(5,4)(2,5)}, {(6,4)(1,5)} 
{(7,5)(2,6)}, {(6,5)(3,6)}, {(7, 6)(4, 7)}, {(0,6)(3,7)} 

The Latin Square to remove above singular fading state 



(figure 21(a) I is obtained by taking transpose of the Latin 

Square shown in figure |20(d)] Then by column shifting all 

, ^ J- -I sin37r/8 

the singular fading states with 7 — 

Next Consider 7 = sin7r/4, f? 
straints are 



is removed. 

smTT/S 

: 0. The singularity con- 



{(0,3)(2,7)}, {(0,5)(6,1)}, {(1,4)(3,0)}, {(1,6)(7,2)} 
{(2,5)(4,1)}, {(3,6)(5,2)}, {(4, 7)(6, 3)}, {(5,0)(7,4)} 



The Latin Square shown in figure |2T(b)] removes this singu- 
lar fade state. Then by shifting columns all the singular fade 
states with 7 = sin7r/4 is removed. By taking transposes all 

the singular fade states with 7 = — — is removed. 

sm 7r/4 

For example, the Latin Square for singular fade state 
7 = sin7r/4, 6* = 7r/2 is obtained by two column shifting 



from Latin Square for 7 = sin7r/4, 6* = 0. The singularity 
constraints for 7 — sin7r/4, = 7r/2 are, 

{(0,1)(2,5)}, {(0,3)(6,7)}, {(1,2)(3,6)}, {(1,4)(7,0)} 
{(2,3)(4,7)}, {(3,4)(5,0)}, {(4, 5)(6, 1)}, {(5,6)(7,2)} 

The Latin Square to remove singular fade state 7 ~ 



sin7r/4, = n /2 is shown in Table 21(c) 



Now consider the singular fade state 7 = —, — , = 

sin 7r/4 

—Tx 12. The singularity constraints for this singular fade state 
are 

{(1,0)(5,2)}, {(3,0)(7,6)}, {(2,1)(6,3)}, {(4,1)(0,7)} 
{(3,2)(7,4)}, {(4,3)(0,5)}, {(5,4)(1,6)}, {(6,5)(2,7)} 

From the Latin Square to remove singular fade state 7 = 
sin7r/4,f? = 7r/2 the Latin Square to remov e singular fade 



1 



2 1 (d) I is obtained by 



state 7 — —, —,9 = — 7r/2 (figure 

sin7r/4 
taking transpose. 

Till now we have removed 40 singular fade states lying on 

the angles 6 — 2m7r/8, where m = 0,1,2- ••7. Now 6 = 

(2m + l)7r/8, where m = 0, 1, 2 • • • 7. is considered. 

The singularity constraints for the singular fade state 7 = 

sin7r/8 

and — 7r/8 are 



sin 47r/8 



(d) 7^ 



sin 37r/8 
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(a) 7 = 


sinyr/S 
sin47r/8' 


e = tt/8 
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(b) 7 = 


sin7r/8 
sin 27r/8' 


6» = 7r/8 
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(c) 7 = 


sin 27r/8 
sin 37r/8 ' 


6» = 7r/8 
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sin7r/8 
sin27r/8 



and 6 — n/8 are 



sin47r/8 

Fig. 22. Latin Squares Corresponding to Different singular fade states 



{(0, 2)(1, 6)}, {(0, 5)(7, 1)}, {(1, 3)(2, 7)}, {(2, 4)(3, 0)} 
{(3, 5)(4, 1)}, {(4, 6)(5, 2)}, {(5, 7)(6, 3)}, {(6, 0)(7, 4)} 

The Latin Square to remove this singular fade state is given in 
figure |22(a)] By shifting the columns of this Latin Square all 
other Latin Squares to remove singular fade states with 7 = 

"^^^ is obtained. The transposes of these Latin Squares 



sin 47r/8 

remove singular fade states with 7 



sin 47r/8 



{(0, 1)(1, 7)}, {(0, 3)(1, 5)}, {(0, 4)(7, 2)}, {(0, 6)(7, 0)} 
{(1, 2)(2, 0)}, {(1, 4)(2, 6)}, {(2, 3)(3, 1)}, {(2, 5)(3, 7)} 
{(3, 4)(4, 2)}, {(3, 6)(4, 0)}, {(4, 5)(5, 3)}, {(4, 7) (5, 1)} 
{(5, 0)(6, 2)}, {(5, 6)(6, 4)}, {(6, 1)(7, 3)}, {(6, 7)(7, 5)} 

The Latin Square to remove this singular fade state is given in 
figure |22(b) By shifting the columns of this Latin Square all 
other Latin Squares to remove singular fade states with 7 = 

— '""^^ is obtained. The transposes of these Latin Squares 
sin27r/8 ^ ^ 

• , -u sin27r/8 
remove singular fade states with 7 = — . 

sin7r/8 

The singularity constraints for the singular fade state 7 = 

sin27r/8 

— and fc* = 7r/8 are 

sin37r/8 ' 

{(0, 2)(2, 7)}, {(0, 3)(2, 6)}, {(0, 4)(6, 1)}, {(0, 5)(6, 0)} 
{(1,3)(3,0)},{(1,4)(3,7)},{(1,5)(7,2)},{(1,6)(7,1)} 
{(2, 4)(4, 1)}, {(2, 5)(4, 0)}, {(3, 5)(5, 2)}, {(3, 6)(5, 1)} 
{(4, 6)(6, 3)}, {(4, 7)(6, 2)}, {(5, 0)(7, 3)}, {(5, 7)(7, 4)} 

The Latin Square to remove this singular fade state is given in 
figure |22(c)] By shifting the columns of this Latin Square all 
other Latin Squares to remove singular fade states with 7 = 

is obtained. The transposes of these Latin Squares 



sin37r/ 



sin 37r/8 

remove singular fade states with 7 = 

® ' sin27r/8 

The singularity constraints for the singular fade state 7 = 

sin37r/8 

and e = V8 are 

{(0, 3)(3, 7)}, {(0, 4)(5, 0)}, {(1, 4)(4, 0)}, {(1, 5)(6, 1)} 
{(2, 5)(5, 1)}, {(2, 6)(7, 2)}, {(3, 6)(6, 2)}, {(4, 7)(7, 3)} 

The Latin Square to remove this singular fade state is given in 
figui'e |22(d)l By shifting the columns of this Latin Square all 
other Latin Squares to remove singular fade states with 7 = 

sin 37r/8 

— ^~"yg obtained. The transposes of these Latin Squares 
remove singular fade states with 7 



L47r/8 



The singularity constraints for the singular fade state 7 



. From above 

sm 37r/8 

discussions it is clear that by obtaining one Latin Square for 
each value of 7 < 1 the other Latin Squares can be obtained 
by column permutation or by taking transpose. The clustering 
for the seven 7 values when end nodes use 8-PSK is explicitly 
given in IV All the clusterings are having 8 clusters only. Thus 
all the singular fade states are removed with 8 symbols only. 
That is, relay can use 8-PSK constellation. 

V. Singularity Removal Constraints FOR M-PSK 
SIGNAL SET 

It was shown in Section III A that the singularity removal 
constraints can be equivalently viewed as a Constrained Partial 
Latin Square (CPLS). In the following subsection, the structure 
of the CPLS for A/-PSK signal set is obtained. 
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Fig. 23. Example showing a 2-plex of order 4 with 8 symbols 
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Fig. 24. Example showing a 2-plex of order 4 
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Fig. 25. CPLS corresponding to the singular fade state 

e Af - e Af sin (^jy j 



A. 77ie Structure of Constrained Partial Latin Square for M- 
PSK Signal Set 

In this subsection, the structure of the CPLS is shown to be 
the one referred to as i-plex which is defined as follows: 

Definition 4: A nx n Partially Filled Latin Square (PFLS) 
is said to be a i-plex of order n with q symbols, if exactly t 
cells are filled in each row and in each column, with each one 
of the q symbols filling exactly nt/q cells in total. A n x n 
PFLS which is a t-plex with n symbols is simply referred to 
as a t-plex of order n. 



For example. Fig. 23 shows a 2-plex of order 4 with 8 
symbols and Fig. 24 shows a 2-plex of order 4. Note that a 
t-plex of order n need not be always completable to a filled 
Latin Square with n symbols. For example, the 2-plex of order 



4 shown in Fig. 24 cannot be completed using 4 symbols, since 
the cell corresponding to the zeroth row and third column 
cannot be filled with any of the four symbols 0, 1, 2, 3. 

The CPLS for the singular fade states which lie on the 
unit circle and the corresponding Latin Squares which remove 
them, were obtained in Section III A. The following lemma 
gives the structure of the CPLS for those singular fade states 
which lie on circles ^r^^-M-r, k ^ I, for the case when both k 
and / are odd or both k and I are even. 



Lemma 6: For a singular fade state 



^' "( M ) „727rm/Af ^JjCrC 



both k and / are odd or both are even, the singularity removal 
constraints are given by. 



{(• 



M [k - I) 
1 2 



k,i — m + 



M 
Y 



jk + l) 
2 



■ k'i — m ■ 



{k - I) 



where < z < M - 1. 

Proof: For {(fci, Zi), (^2, ^2)} to be a singularity con- 
straint corresponding to the singular fade state " ) , 



. 7r(fci -A:2) 



sin( 



M 

n(l2-h) 
M 



sin 




gj27rm/M 


sin 







(8) 



or. 



Equating the amplitude and phase terms in (|8]l, we get, 

ki + k2 = h + h, ki — k2 = k,l2 — h = I 

ki + k2 = h + h + 2m, fci — fc2 = k, I2 — h = M — I. 



Note that by Lemma [T] separately equating the numerator 
and denominator of the amplitude on both the sides of ([8]) 
is allowed. 

Assuming fci = i, the two solutions given in the statement 
of the lemma are obtained. ■ 
For example, the CPLS corresponding to the singular fade 



state I ? , is shown in Fig 



25 



TITe structure of the CPLS for 



sin(f) 

The following lemma gives 
those singular fade states which lie on circles jl^ ,k ^ I, 
for the case when only one among k and I is odd. 



Lemma 7: For a singular fade state 



where only one among k and I is odd, the singularity removal 
constraints are given by. 



{( 
{( 



M 
2 

(fc- 



(fe + 1 - 



-0 



~ k,i 
k,i — m - 



(k 



M 

1-0 



(fc + 1 + 



)} 



)}^ 



where < i < M - 1. 

Proof: The proof is similar to the proof of Lemma [6] and 
is omitted. ■ 

The structure of the CPLS for A/-PSK signal set, is shown 
to be a 4-plex of order M with 2M symbols or a 2-plex of 
order Af, as stated in the following lemma. 

Lemma 8: The CPLS obtained after filling in the singularity 
removal constraints corresponding to the singular fade state 



which lies on the circle with radius 



fc 7^ Z, forms a 

sin(^)' ' 

4-plex of order M with 2M symbols if fc 7^ M/2,/ 7^ M/2 
and forms a 2-plex of order M otherwise. 

Proof: Case (i) k ^ M/2, I ^ M/2: From Lemma |6] 
and Lemma [7] each value of i, < i < M — 1, gives rise 
to two distinct constraints, which fills two cells in the z-th 
row. Similarly, the two constraints corresponding to i' = i + k 
also fills two cells in the i-th row. Hence, in each row, 4 cells 
are filled. By a similar argument, 4 cells are filled in each 
column. Since the same symbol fills in two of the cells, a 
total of 4M/2 — 2M symbols are used in the partial filling 
of Latin Square. As a result, a 4-plex with 2M symbols is 
obtained. 

Case (ii) k = M/2 or I ^ M/2: When k = M/2 or I = M/2, 
from Lemma |6] and Lemma |7] for each value of i, < « < 
M — 1, it can be verified that the two singularity constraints 
are not distinct. Only 2 cells are filled in each row and each 
column. As a result, a 2-plex is obtained. ■ 
For example, the PFLS corresponding to the singular fade 



shown in Fig. 



25 



is a 4-plex of order 8 with 16 



state . 
symbols. 

If the same symbol is filled in those cells which belong to 
two different singularity constraints, the singularity removal 
constraints are said to be combined. It is shown in the 
following lemma that by combining the singularity constraints, 
the CPLS which is a 4-plex of order M with 2M symbols can 
be transformed into a 4-plex of order M. 

Lemma 9: For k ^ M/2, 1 ^ M/2, the CPLS corre- 
sponding to the singular fade state which lies on the circle 
'^!"^ \ ,k^l, obtained after filling in the singularity con- 
straints, can be transformed in to a 4-plex by appropriately 
combining the singularity constraints. 

Proof: Since the elements of both the sets — k,i + 
M/2,i-k + M/2} and {i-M/2~ik~l)/2,i + M/2-{k + 
l)/2, i - {k - l)/2,i - {k + l)/2} are distinct, the following 
constraints in Lemma |6] can be combined: 
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i,i 
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i H , i - 

2 
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(k-l) 



k,i + 



i — k - 



M (k + I) 
"2 2 

(k + l) 



M 
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2 ' 



Similarly, the following constraints can also be combined: 



(k + I) 



M M 

i-\ ,«H 

2 2 



i — k,i 
{k + l) 



(k-l) 



M M 

fcH ,i+ — 

2 2 



(k-l) 



In this way, a 4-plex is obtained from the PFLS for the 
case when k,l ^ ^ and both k and I are even or both are 
odd. Following, exactly the same procedure, the lemma can 
be proved for the case when one among k and I is odd and 
the other is even. ■ 
For example, the CPLS corresponding to the singular fade 



(^) 



state "'"\^/ , shown in Fig. 25 
plex of order 8, as shown in 



can be transformed into a 4- 
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Fig. 26. 4-plex of order 8 corresponding to tlie singular fade state 



Fig. 26 



The completability of a <-plex of order M with M symbols 
is discussed in the following subsection. 



B. Completability of a t-plex of order M with M symbols 

The completion of the 2-plexes and the 4-plexes of order 
M, obtained from the singularity removal constraints, using 
exactly M symbols is advantageous for the reason outlined 
as follows: It was shown in Q that for the case when QPSK 
signal set is used during the MA phase, there exists certain 
values of fade state for which a clustering with cardinality 5 
needs to be used to maximize the minimum cluster distance. 
While the use of clusterings with cardinality 5, reduces the 
impact of multiple access interference, it adversely impacts 
the performance during the BC phase. The 2-plexes and 4- 
plexes can be completed using M symbols, means that the 
optimizing the performance during the MA phase does not 
come at the cost of a degraded performance during the BC 
phase. 

Whether a partially filled M x M Latin Square is com- 
putable using M symbols or not is an open problem. Burton 
Q conjectured the following regarding the completability of 
the partially filled Latin Square which forms a <-plex of order 
M. 

Conjecture 1: A t-plex of a Af x M Latin Square is 
completable using M symbols, for t < M/4. 

By explicitly providing all the Latin Squares, for 8-PSK 
signal set, it was shown in the previous section that the CPLS 
corresponding to all the singular fade states are completable 
using 8 symbols. If Conjecture [T] were indeed true, the 4-plexes 
and the 2-plexes obtained from the CPLS corresponding to the 
different singular fade states, can be completed using M sym- 
bols, for M any power of 2 greater than or equal to 16. Even 
though, the problem of completability of a i-plex in general 
remains unsolved, in the next section, it is shown that the 
structure of the singularity constraints allows the construction 
of explicit Latin Squares, for some singular fade states. In 
other words, by providing some explicit constructions, it is 
shown that the 4-plexes of order M corresponding to some 
of the singular fade states are always completable using M 
symbols. 



VI. SOME SPECIAL CONSTRUCTIONS OF LATIN 
SQUARES 

As discussed in section II, the singular fade states lie on 



circles with radii "'" / t',\ ,i < k,l < M/2. It was shown in 



Section II that all the Latin Squares which remove the singular 
fade states on the unit circle (i.e., the case when k = I) can be 
obtained from the bit-wise XOR map, by appropriate column 
permutation. In Subsection A of this section, an explicit 
construction algorithm is provided to obtain Latin Squares for 
the case when both k and I are odd. In Subsection B of this 
section, a construction procedure to obtain Latin Squares for 
the case when both k and I are even is provided. Note that the 
cells in the Latin Squares obtained are filled using exactly M 
symbols. 

It was observed in |4| that for 4 PSK signal set there exists 
clusterings which remove more than one singular fade state. 
The Latin Squares obtained from the construction algorithm 
given in Subsection A remove multiple singular fade states. 

The construction procedure given in Subsection B involves 
obtaining Latin squares for Af-PSK signal set from the ones 
already obtained for M /2-PSK signal set. 

A. Construction of Latin Squares which remove Multiple Sin- 
gular Fade States 

Let £^ I denote the Latin Square which removes multiple 



singular fade states on the circle 



,1 < k,l < M/2, 



both k and I are odd, which is obtained as follows: 
1: Start with an empty Latin Square. 

2: The first row of the M x M Latin Square is filled with 

numbers to M — 1 in the same order. 
3: for < s < A/ - 1 do 
4: for 1 < i < M do 
5: if s is even then 

6: The entry {kt, s + It) is filled in the Latin Square 

is filled to be s. 
7: else 

8: The entry {kt, s — It) is filled in the Latin Square 

is filled to be s (Note that all the arithmetic 
operations are modulo M.) 
9: end if 
10: end for 
11: end for 

Since both k and I are odd, the algorithm given above is always 
guaranteed to result in a completely filled Latin Square. 

Let I denote the Latin Square obtained by replacing the 
word even by the word odd in Step 5 of the algorithm. 

For example consider the case when k ^ 3,1 — 1. We start 



with only the first row filled as shown in Fig. 27(a) The cells 



in which O's appear are filled next, by traversing the Latin 
Square, fc = 3 steps forward and I = 1 step to the right, as 
depicted in Fig. |27(b)l The cells in which I's appear are filled 
next, by traversing the Latin Square, fc = 3 steps forward and 
I = 1 step to the left, as depicted in Fig. 27(c) In a similar 



way, the remaining cells in the Latin Square can also be filled 
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(b) 
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Fig. 27. Construction of the Latin Square £| 



Fig. 28. The Latin Square £g j 



to obtain the Latin Square L\ shown in Fig. 27(d) Similarly, 
the Latin Square £3 can be obtained as shown in Fig. 28 

The following Theorem states that multiple singular fade 
states on multiple circles are removed by the Latin Squares 
l^l.i and Cl ^. 

Theorem 1: One of two the Latin Squares £| ; and ; 
(both k and I odd), removes any singular fade state of the form 

sin '^l^ ym1-K 

— j^f^eT?^, where < m < A/ — 1 and (k\.li) belongs to 

sin ^ 

the set, {{nk,nl)\\ < n < ^ - l,n odd}. 

Proof: The singularity constraints for the singular fade 
state f'i e^T"^ are given by Lemma 6 

Assuming i' = i — k. The singularity constraints given in 
Lemma |6] can be rewritten as. 



2 2 / V 2 2 



(9) 



(fe + 



i' + k,i' + m + 



(k - I] 



(10) 



where < i' < M - 1. 

From (|9]l, the second entry in the the singularity constraint 
can be obtained from the first entry by traversing k rows 
downwards and I columns to the left. Similarly, from ( fTO] ), the 
second entry in the the singularity constraint can be obtained 
from the first entry by traversing k rows downwards and I 
columns to the right. Furthermore, for a given i', the column 
indices of the first cells given in the singularity constraint given 
in (|9]) and ( [TOj l differ by an odd number. Hence, one of the 
two Latin squares £^ ^ and ; removes the singular fade state 



, depending on whether m 



(fc+0 



is odd or even. 



From, the construction procedure described in the beginning 
of the sub-section, it can be verified easily that C^. ; is the 
same as £^j, and £^ ^ is the same as £° j, for all odd n. 
This completes the proof. ■ 
It follows from Theorem |T] that a single Latin Square Cf, ; 
removes half of the singular fade states which lie on M/4 
circles, i.e., a total of ^ singular fade states. Two Latin 
Squares C 



k.l 



and 



remove all the fade states which lie 



on M/4 circles. Also, the Latin square £^ ^ (£JJ j) is the same 

iKk.ni\ where3<n<A//2-l, 



as the Latin Squares C'^^. 
n odd. 

For example, the Latin Squares £3 ^ and £3 
respectively in Fig. 27(d) and Fig 



shown 

remove all the 16 
singular fade states which lie on the 2 circles with radii 



28 



sm( - 



, where (fci,/i) e {(1, 3), (3, 1)}. 



Note 2: The Latin Square £^ ^ is the same as the one 
obtained by a single column permutation of £^ with the 
symbols alone permuted. Since the clustering is unchanged 
by symbol permutation, it is enough if we specify either Cj. ; 



Out of the 



4 



M 
2 



States lie. 



1 circles on which the singular fade 
44 ^ circles have radii of the form ^ifl} , 

16 4 sin(ij)' 

where both k and I are odd. From Theorem 1, the Latin 
Squares which remove all the singular fade states which lie 



on such circles can be obtained from 



= ^ - 1 Latin 



Squares, which are constructed using the algorithm given in 
the beginning of this subsection. 

B. Construction of Latin Squares for M PSK signal set from 
that of M/2 PSK signal set 

In the previous subsection, an explicit construction of Latin 
Squares which remove the singular fade states on the circle 
^™ ^ , both k and I odd, was provided. In this subsection, 

sin -jj 

it is shown that for M PSK signal set, Latin Squares which 
remove the singular fade states on the circle ^™ , both k 
and I even, k ^ l,k ^ M /2 and I ^ M/2, can be constructed 
from the Latin Squares obtained for M/2 PSK signal set. 
Let C{M, k, I, (j)) denote an M x M Latin Square which 

removes the singular fade state fL I e^'^. 

For the M x M Latin Square L{M,k,l,(j)), let 
C°''{M, k, I, (j)) denote the M /2 x M /2 Latin Square obtained 
by taking only the odd numbered rows and even numbered 
columns of £(M, fc, I, 0). In a similar way, the the M/2 x M/2 
Latin Squares £'=^(M, fc, Z, </)), C°°[M, fc, Z, </>), £'=°(M, k, Z, 0) 
can be defined. 

Definition 5: The quadruplicate of an M x M 
Latin Square C(AI, k,l, (f>) is the ordered set 
consisting of the four M/2 x M/2 Latin Squares 
{C^^M, k, I, 0), £°°(M, k, I, 0), £'=°(M, k, I, (/>),£°'=(M, k, I, (j))}. 

The Latin Square is C{M, k,l, (f>) uniquely determined by 
its quadruplicate. 

For the Latin Square C{M, k, I, (j)), let Cc{M, k, I, (p) denote 
the Latin Square obtained by adding integer c to all the cells 
of the Latin Square. 

The following theorem gives the construction of the Latin 
Square £(M, fc, /, (j)), for the case when both k and I are even, 
k^l,k^ M/2, I ^ M/2 and mod (| + |, 2) = 0. 

Theorem 2: The elements of the quadruplicate of 
£(M, fc, I, cj)), if mod (I + |, 2) = 0, ai-e given by. 



C^iM, k, /, (t>) = C°°{M, k, I, </>} = CiM/2, k/2, 1/2, </>), 
£=°(M, k, I, (p) = C°''{M, k, I, <)>) = Cm/2{M/2, k/2, 1/2, 



Proof: From Lemma |6j the singularity constraints are 
given by, 



M (k-l) 
~2 2 



M (k + l)\ , 

(11) 



(12) 



where < i < M - 1. 

The singularity constraints given above can be further split 
into cases depending on whether i is even or odd. 
Case 1: i is even 

Assume i' = i/2, the singularity constraints in (H) and 
(fT2ll can be rewritten as. 



2j',2i' 



M/2 2{k/2- 



2i' - 2-,2j' - 



(k/2 + l/2) 



^M/2 Jk/2 + 
2 2 



1/2) 



2i' 



k I 
2 



(k/2-l/2) 



)}■ 



(13) 



(14) 



where < i' < M/2 - 1. 

The singularity constraints given in ( [T3] l and <^4]), are the 



same as that of the singular fade state 



M/2 I 



with the 



row and column indices doubled. Hence, C{M, 2k, 21, (jyy^ = 
C{M/2,k,l,(j)). 
Case 2: i is odd 

Using a similar argument as that of case 1, it can be 
shown that the singularity constraints for this case are the 

sin( -f-TjTT ) 

same as that of the singular fade state with the 



. / Ttl/2 s 

>5in(M72) 



row and column indices doubled and shifted by one place. 
Hence, C°°{M,k,l,(i)) = C{M/2,k/2,l/2,(j)). Note that 
by taking both C"'{M,k,l,(t>) and C°°{M,k,l,(j)) to be 
L[M/2, k/2,l/2, <j)), the conditions required for £(M, k, I, 4>) 
to be a Latin Square are not violated. 

Since C^<'{M,k,l,4>) = C°°{M,k,l,4>) = C{M/2,k/2,l/2,4>), 
half of the cells in the Latin Square C{M,k,l,<f) are filled 
with numbers from to M/2 - 1. Hence, C°''{M, k, I, 0) and 
^"{M, k, I, (f)) can be taken to be Cm/2{M/2, k/2,l/2, cj)). 

■ 

For example, the Latin Square £(16,6,2,0) can be ob- 
tained from £(8,3,1,0) (which is shown in Fig. 28 1 as 
follows: £'=•=(16,6,2,0) and £5°°(16, 6 2, 0) ai-e chosen 
to be £(8,3,1,0) as shown in Fig. 29(a) The Latin 



Squares £°'=(16, 6, 2, 0) and £'=°(16, 6, 2, 0) are chosen to be 
£8(8,3,1,0) to obtain the Latin Square £(16,6,2,0) shown 
inFig.|29(b)l 

The following theorem gives the construction of the Latin 
Square £(M, k, I, </)), for the case when both k and I are even, 
k^l,k^ M/2, I ^ M/2 and mod (| + ^, 2) 1. 

Theorem 3: The elements of the quadruplicate of 

£(M, k, I, 4)), if mod (| + |, 2) = 1, are given by. 
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Fig. 29. Construction of the Latin Square £(16, 6, 2, 0) 



£°=(M, k, I, <f>) = C{M/2, k/2,l/2, 0), 

£'=°(M, k, I, 0) = C{M/2, k/2,l/2, 4> - 2-k/M), 

C°°{M, k, I, (/)) = C^^iM, k, I, (f)) = Z;j,,^/2(M/2, fe/2, Z/2, ij>). 

Proof: The proof is similar to Theorem |2] and is omitted. 

■ 

In this subsection, construction procedure was provided to 
obtain the Latin Squares which remove the singular fade states 



which lie on circles with radii 



where both k and I 



are even, k ^ l,k ^ M/2 and / ^ M/2. The total number 



of such circles is (4^ — 1 



2. 



( K — ^\ — Ml _ 3M 
It is enough if we obtain one Latin Square per circle for 
those circles which lie inside the unit circle, since the Latin 
Square corresponding to all other singular fade states can be 
obtained by appropriately performing column permutation and 
taking transpose. Hence, to remove all the singular fade states 



sin(^ 



where both k and I are even. 



on circles with radii . , , . 

sin(lif ) 

,k 7^ M/2 and I ^ M/2, it is enough if we obtain 



k ^ I, 



2 Latin Squares from the construction procedure 

described in this subsection. 

The Latin Squares constructed so far (the bit-wise XOR map 
and the ones obtained in Subsection A and B of this section). 



which are 



'M 



16 



3M 
4 



M2 Af 

+1 

32 8 



76- 



in number, are sufficient to obtain all the Latin Squares which 
remove the singular fade states which lie on 



1 + 



M2 

16 



M 

T 



3M 



/M2 

V 16 4 



M2 



M + 3 



circles. The Latin Squares which remove the singular fade 
states which lie on the other circles inside the unit circle can 
be found using computer search. It can be verified that it is 

^ Latin Squares in total, from 



sufficient to obtain ^^||- 



which all the other Latin Squares can be obtained, for M > 8. 

VII. Bidirectional Relaying With Unequal 
Transmission Rates and Latin Rectangles 

In this section we consider the scenario in which both 
the end nodes use PSK constellations of different sizes. This 
results in different transmission rates for both the users. 
Assume node A and node B use PSK constellations Sa and 
5b of size M{= 2'^) and N{— 2^) respectively. The relay 
should use a constellation of size at least niax{AI,N). The 
exclusive law for this case is given by. 



M~''^{xa,xb) ^ M~''^{x'^,xb), where xa x\ , xb £ Sb, 
M~''^(xa,xb) ^ M~''^{xA,x'g), where xb x'g , xa & Sa- 



(15) 



Definition 6: A Latin rectangle L on the symbols from the 
set Zt = {0, 1, • • • , t — 1} is an M X array, in which each 
cell contains one symbol and each symbol occurs at most once 
in each row and column. 

The relay clusterings which satisfy the exclusive law form 
Latin rectangles with M rows and N columns. As before, 
the study of clusterings which satisfy the exclusive law can 
be equivalently carried out as the study of Latin rectangles 
with appropriate parameters. Without loss of generality, it is 
assumed that M > N. It is shown that the set of singular fade 
states for this case is a subset of the set of singular fade states 
for the case when both A and B use 7\/-PSK signal set. Also, 
it is shown that the M x N Latin Rectangle for this case can 
be obtained from M x M Latin Square obtained for the case 
when both A and B use A/-PSK signal set. 

Lemma 10: The set of singular fade states for the case when 
Af-PSK signal set is used at A and A^-PSK signal set is used 
at B, both M and N powers of 2 and M > N,\s a. subset of 
the set of singular fade states for the case when both A and 
B use M-PSK signal set. 

Proof: Let Sm and Sn denote the symmetric A/ -PSK 
and A^-PSK signal sets used at A and B respectively. 

Assume M = 2^ and = 2^, A > 

Note that the signal set Sjsi is a subset of the signal set Sm- 
The singular fade states for this case are of the form. 



x{k) - x{k') 

' y{i)-y(i') 



x{k) ^ x{k') e 5m, 2/(0 7^ y{l') e Sn (16) 



The singular fade states given by ([T6|, are also singular fade 
states for the case when both A and B use i\/-PSK signal set, 
since y{l) and y{l') are points in the il/-PSK signal set as 
well. This completes the proof. ■ 

The M X N system can be interpreted as a M x M system 
with node B not transmitting some symbols. The correspond- 
ing Latin Rectangle is obtained from the Latin Square for 
M X M system by removing the columns corresponding to 
the unused symbols. The following example illustrates this. 

Example 5: Consider the case when 8-PSK signal set is 
used at A and 4-PSK signal set is used at B. For 7 = sin(7r/4) 
and 6* = 0, by removing the columns 1, 3, 5 and 7 (equivalently 
viewed as B transmitting only the 4 symbols 0, 2, 4 and 6 in 
the 8-PSK signal set) from the Latin Square given in Fig. 
|21(b)| the Latin Rectangle shown in Fig. |30]is obtained. 
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Fig. 30. Latin Rectangle for 7 = sin(7r/4), 6 = 



VIII. DISCUSSION 

In this paper, for the design of modulation schemes for 
the physical layer network-coded two way relaying scenario 
with the protocol which employs two phases: Multiple access 
(MA) Phase and Broadcast (BC) phase, we identified a relation 
between the required exclusive laws satisfying clusterings and 
Latin Squares. This relation is used to get all the maps to 
be used at the relay efficiently. Further we illustrated the 
results presented for the case, where both the end nodes use 
QPSK constellation as well as 8-PSK constellations. Here we 
concentrated only on singular fade states and the clusterings 
to remove that with only the minimum cluster distance under 
consideration. We are not considering the entire distance 
profile as done in |4|. Our work eliminate the singular fade 
states effectively and these clusterings can be used in other 
regions in the complex plane of (7, 9), as shown in pO) . 
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